Kelvin waves of quantized vortex lines in trapped Bose-Einstein condensates 
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We have theoretically investigated Kelvin waves of quantized vortex lines in trapped Bose-Einstein 
condensates. Counter-rotating perturbation induces an elliptical instability to the initially straight 
vortex line, driven by a parametric resonance between a quadrupole mode and a pair of Kelvin modes 
of opposite momenta. Subsequently Kelvin waves rapidly decay to longer wavelengths emitting 
sound waves in the process. We present a modified Kelvin wave dispersion relation for trapped 
superfluids and propose a simple method to excite Kelvin waves of specific wave number. 
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In classical hydrodynamics a vortex line is the trajec- 
tory obtained by following the local vorticity vector, de- 
fined as the curl of the velocity field of the fluid [1] . Wing 
tip vortices of aircrafts, tornadoes, whirlpools in a flow- 
ing water, and cosmic strings provide intuitive mental 
pictures of vortices. They support a branch of chiral nor- 
mal modes in which the perturbation propagates along 
the vortex line and rotates about its unperturbed posi- 
tion therefore distorting the vortex into a helical shape. 
In 1880, Thomson (Lord Kelvin) derived the dispersion 
relation for such excitation modes, which are now known 
as Kelvin waves 0. These helically vibrating normal 
modes are objects of fundamental importance in classi- 
cal hydrodynamics. 

Similar structures also exist in superfluids in which vor- 
ticity is quantized in units of h/m, where h is Planck's 
constant and m is the mass of the superfluid particle 
[H S EL HI, @i 0] ■ Kelvin modes or kelvons of singly quan- 
tized rectilinear superfluid vortices have no radial nodes 
and in axisymmetric systems they are characterized by 
their angular momentum quantum number I = — 1 and 
momentum Kk z along the vortex axis. First experimental 
evidence of Kelvin waves in superfluids was obtained by 
Hall and they are thought to play an essential role in the 
formation and decay of superfluid turbulence [8| • In par- 
ticular, Kelvin wave and Kolmogorov-type cascades are 
predicted to mediate the energy transfer across length 
scales in turbulent superfluid systems 0, QjJ II , H, 3 • 
Bretin et al. were able to excite and observe Kelvin waves 
in a single quantum vortex in Bose-Einstein condensates 
[T3 |. Their experimental kelvon excitation mechanism 
was theoretically clarified by Mizushima et al. 
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The ability to address individual vortex lines in these 
systems opens the possibility to study Kelvin waves and 
their relation to elementary processes governing super- 
fluid turbulence. 

In this Letter we study the microscopic dynamics of 
Kelvin waves of quantized vortex lines in Bose-Einstein 
condensates by employing direct simulations of the three- 
dimensional time-dependent Gross-Pitaevskii equation 
combined with an analysis of the collective Bogoliubov 
excitations of the unperturbed state. We treat the dy- 



FIG. 1: (Color online) Kelvin waves in a quantum mechan- 
ical vortex. The initially straight vortex (a) undergoes an 
elliptical instability and the vortex line shape changes from 
plane-sinusoidal (b) to helical (c). Subsequently the primary 
kelvons decay into kelvons of longer wavelength (d) emitting 
sound waves seen as ripples at the condensate surface. In the 
laboratory frame the sense of the vortex motion changes sign 
between frames (c) and (d). The frames are for times ms (a) 
103 ms (b) 129 ms (c), and 350 ms (d). The bottom end-cap 
of each frame has been cut-off for the aid of visualization. 



namics of Kelvin waves fully quantum mechanically in 
contrast to previously employed phenomenological mod- 
els such as in Refs.[17l. |l8j. We find counter-rotating 
quadrupole surface mode (surfon) excitations to decay 
via Beliaev mechanism into a pair of kelvons of oppo- 
site momenta in perfect agreement with Refs. 
Subsequently the primary kelvons are found to decay into 
kelvons of longer wavelength by emitting sound waves in 
the scattering process. This mechanism may help to ex- 
plain microscopically how superfluid turbulence returns 
to laminar flow even at zero temperature. We introduce 
a modification to the semi-classical kelvon dispersion re- 
lation extending its validity to trapped superfluids. Fi- 
nally, we propose an experimentally feasible method to 
excite kelvons of particular wave number. 

Following the experimental procedure of Bretin et al. 
[H | we consider a Bose-Einstein condensate composed of 
1.3 x 10 5 87 Rb particles trapped in a harmonic potential 
whose transverse and axial frequencies are u>±_ = 2ir x 98.5 
Hz and uj z = 2tt x 11.8 Hz, respectively. We model the 
physics of such system by numerically solving the pure 
(no additional damping terms included) time-dependent 
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Gross-Pitaevskii equation on a parallel computer us- 
ing the method described in Refs. We com- 
plement the dynamical simulations by directly solving 
the Bogoliubov-de Gennes eigenvalue equations obtain- 
ing the full Kelvin wave dispersion relation [3, [l6[ . Our 
initial condition is a near-axisymmetric single-quantum 
vortex which possess one unit of angular momentum per 
particle. Such state, shown in FigQJa), is a local energy 
minimum in the non-rotating trap rather than an abso- 
lute ground state which would have zero vorticity. 

At t = we suddenly switch on an elliptical rotating 
perturbation 

V pcrt {r,t) = ^emu; 2 _[{x 2 - y 2 ) cos(2ilt) + 2xy sin(2fii)], 

(1) 

where e = 0.025, and m is the mass of the atom. The fre- 
quency f2 = — 0.6wj_ is resonant with the counter-rotating 
(with respect to the superfluid flow) £ = — 2 quadrupole 
surfon of the unperturbed condensate This reso- 

nant external drive transfers population from the initial 
vortex state to the counter-rotating quadrupole surfon 
state thus changing the total orbital angular momentum 
of the system. When the perturbation is switched off 
at t = 35ms, the ^-component of the angular momen- 
tum has been reduced from L z /N = lft to 0.47ft. From 
there on L z is a constant of motion since the bare trap 
potential is cylindrically symmetric. In contrast, if the 
system is continuously rotated, the growing surface insta- 
bility leads to (a nti) vortex nucleation as in nonrotating 
systems [2l], 22, 23, 24, [25| and eventually Lz changes 
its sign. If £1 is held constant the angular momentum 
transfer is limited since the rotating perturbation falls 
out of resonance when the quadrupole surfon becomes 
frequency shifted due to the changing angular momen- 
tum. 

The vortex line stays straight until t — 64 ms at 
which stage the quadrupolar condensate density modula- 
tion in the plane perpendicular to the vortex axis induces 
squeezing in the vicinity of the vortex core. The super- 
flow then undergoes an elliptical instability, the vortex 
line length stretches and its shape deforms from straight 
to sinusoidal, oscillating in the plane of the elliptical de- 
formation as illustrated in FigQJb). The elliptical in- 
stability explains how the two-dimensional perturbation 
yields three-dimensional flow instability, and it is driven 
by a parametric resonance between the I = — 2 surfon 
and a pair of I = —1 kelvons with opposite axial mo- 
menta. The elliptical instability mechanism has been 
reviewed by Kerswell in the context of classical hydro- 
dynamics [26j . As the instability grows the shape of the 
vortex changes from plane-sinusoidal to helical as shown 
in FigQTc). Both the vortex core radius and the am- 
plitude of these Kelvin waves are larger near the ends 
of the condensate, explained by the lower local parti- 
cle density in those regions. These Kelvin waves decay 
rapidly via phonon emission into kelvons of longer wave- 
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FIG. 2: (Color online) Kelvin wave excitation and decay (a) 
and damping of quadrupole modes (b) as functions of time. 
Frame (a) is obtained by Fourier transforming the azimuthal 
angle of the vortex position vector as a function of z. The 
quadrupole moments in (b) are computed for £ = —2 (solid 
line) and £ — +2 (dashed line) perturbations, respectively. 
The envelope functions are exponentially vanishing with the 
decay constant F = 16s -1 . 



length. The phonons are visible as ripples on the con- 
densate surface in Fig[TJ After the primary kelvons have 
decayed the system is left in a bent vortex state, FigQJd), 
and the final configuration is determined by the amount 
of angular momentum left in the system after the ex- 
ternal perturbation is switched off. If, instead, the co- 
rotating quadrupole surfon is resonantly populated using 
il = 0.8u>± in Eq. (fTJ), the vortex remains straight for 
all times and the £ = +2 surfon population is undamped 
as is verified by our simulations. In this case the con- 
served angular momentum of the system after the 35 ms 
excitation is L z /N — 1.51ft and the condensate density 
oscillates with the characteristic quadrupole frequency 
/ fc+2 = 162Hz. 

To further quantify the kelvon excitation and decay 
process, we have located the planar position (x v ,y v ) 
of the vortex phase singularity as functions of z and 
t. Fourier transform T of the complex signal Q(z) — 
a,rctaii(y v /x v ) + i a,rcta,n(x v / y v ) reveals clear peaks at 
wave vectors corresponding to the kelvon excitations 
present in the system. Figure Ufa) shows |J r [9(z)]| 2 
as a function of time for the £ = — 2 initial perturba- 
tion. There emerges a strong signal at |k z | = 0.9/im" 1 
which shifts rapidly around 150 ms to |k' 2 | = 0.2/im -1 . 
Figure [2jb) shows the normalized quadrupole moment 
Q(t) = J xyn(r) dr, where n(r) is the condensate den- 
sity, of the system as function of time for the I = — 2 
(solid line) and £ = +2 (dashed line) perturbations. The 
signal for counter-rotated case is strongly damped due to 
the resonant surfon-kelvon coupling while the co-rotated 
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FIG. 3: (Color online) Kelvin wave dispersion relation show- 
ing Bogoliubov excitation frequencies of the I — phonon 
(circles) and £ = — 1 kelvon (triangles) modes as functions 
of the axial momentum. The square at u) = 1.21^x denotes 
the I — — 2 surfon and the two diamonds near the kelvon 
dispersion are the initial and final kelvon modes whose wave 
vectors and energies are obtained from the dynamical simu- 
lations. The solid line and curve respectively correspond to 
the analytical phonon and kelvon dispersion relations in the 
text. The dashed curve demonstrates the failure of the origi- 
nal kelvon dispersion relation. 



surfons remain essentially undamped. The exponentially 
decaying envelope oc e~ rt , with T = 16s" 1 in Fig[^b) il- 
lustrates the damping of the quadrupole oscillation with 
the carrier frequency ft=— 2 = 124 Hz. Landau damp- 
ing due to the 30% thermal fraction in the experiment 
is likely to explain the faster observed damping rate of 
57 ± 10s" 1 Q. 

To complete the picture, we have computed the long- 
wavelength Kelvin wave dispersion relation by explicitly 
solving the Bogoliubov-de Gennes equations for the ini- 
tial state corresponding to FigOJa) [3, lj| 



27]. Figure 



[3] shows the excitation frequencies for odd (open mark- 
ers) and even (solid markers) z-parity modes as func- 
tions of axial momentum for I = — 1 kelvons (trian- 
gles) and I = phonons (circles). The straight line is 
plotted for lu p = c z \k z \, where the chemical potential fi 
in the speed of sound c = \Jb\ijm is renormalized by 
the factor b = 1/3 due to the inhomogeneous density. 
The square at u = 1.21a;_i_ denotes the counter-rotating 
I = — 2 quadrupole excitation. This surfon state is ini- 
tially populated by the resonant (Q = — 0.6uj±) pertur- 
bation. The two lowest energy kelvon modes are the de- 
generate "U" (even parity) and "N" (odd parity) shape 
bending modes [28} while the second lowest modes corre- 
spond to the even parity "£", seen in FigQJd), and odd 
parity "Z" shapes (see also FigHJ. Dynamically these 
in-plane bending modes preserve their shape and they 
precess about the z-axis of the trap. The presence of the 



FIG. 4: (Color online) Bogoliubov kelvon modes for lowest 
even (left) and odd (right) z-parity modes corresponding to 
those in Fig|3] The numbers on top and bottom of each frame 
indicate the mode frequency and momentum, respectively. 
The color indicates amplitude and sign of the displacement 
from the equilibrium position 



bending modes in the spectrum manifests the inhomo- 
geneous density due to the trapping. Notice also that 
the axial momentum is not a good quantum number due 
to the axial trap potential. Nevertheless, we may as- 
sign a well-defined momentum fck z for the Bogoliubov 
modes [iH] . The two diamonds near the kelvon dispersion 
are the initial and final kelvon modes, shown respectively 
in FigQJb) and Figfljd), whose wave vectors and ener- 
gies (rotation frequencies) are extracted from the time- 
dependent simulations. The sign of the kelvon frequency 
determines the sense of rotation of the vortex in the lab- 
oratory frame. In FigU]we have plotted the Bogoliubov 
itfc(x, z) modes for a few of the lowest collective even (left) 
and odd (right) z-parity kelvon excitations 15, lq . The 
numbers on top and bottom of each frame indicate the 
excitation frequency and momentum, respectively, of the 
mode and correspond to those in FigJ3] All kelvons are 
highly localized within the vortex core. 

The Kelvin wave dispersion relation involving Bessel 
functions is derived for an infinitely long classical vortex 
and its long wavelength limit is often applied to super- 
fluid vortices by simply replacing the classical circulation 
by its quantum mechanical counterpart Guided by 
our numerical experiments for N in the range (1 — 50) x 
10 4 , we propose a modification to such semi-classical dis- 
persion relation extending its validity to the finite and 
inhomogeneous systems: 



w(fc 



LO 



2m 



log 



1 



c k 2 |<l (2) 



where ujq is the frequency of the kelvon with the small- 
est axial momentum hko and the system specific vortex 
core parameter r c = 0.13/xm is only weakly dependent on 
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the density. As shown in Fig. [3l Eq. ([2]) (solid curve) 
is in an excellent agreement with the Bogoliubov spec- 
trum (except for the pure bending modes). We have also 
plotted Eq. (T5|) using lu = k = and r c = £ = 0.3/im 
(dashed curve) showing the failure of the original kelvon 
dispersion relation. We emphasize that the usual den- 
sity dependent healing length £ provides poor estimate 
for the constant r c for all values of N studied here. 

Here we propose a simple method to excite kelvons of 
specific wave number on demand. We may still make 
use of the parametric surfon-kelvon resonance by simply 
shifting the whole kelvon dispersion curve with respect 
to the quadrupole frequency. This is achieved by adding 
a Gaussian pinning potential 



Vp in (a;, y, z) = _ / _ ^ exp 



-2[x 2 



y 



(3) 



where <j{z) = cros/l + (z/zr) 2 and we have chosen = 
2/im and a Rayleigh range of zr = 20/xm. The vor- 
tex pinning effectively modifies the local potential at the 
vortex core [2l], |2^] . The strength of the applied poten- 
tial Vq may be used to control the wave vector of the 
kelvon resonant with the quadrupole surfon. Since the 
pinning potential is localized in the centre of the trap 
it does not radically affect the resonant quadrupole sur- 
fon frequency itself. To demonstrate the feasibility of 
this method we have performed simulations using differ- 
ent values for Vq in Eq. ([3]) the only difference to the 
previously described excitation method being the addi- 
tion of the pinning potential at t — 0. Setting Vq ={8.0, 
6.0, 4.0, 2.0, -2.0, -4.0}fio;_L results in an initial kelvon 
excitation with |k 2 | ={0.2, 0.4, 0.7, 0.7, 0.9, 1.0} fmr 1 , 
respectively. Experimentally, this method could be real- 
ized using tightly focused optical potentials allowing full 
control over the kelvon excitation process. 

In conclusion, we have theoretically studied the mi- 
croscopic excitation and decay mechanism of superfluid 
Kelvin waves in trapped Bose-Einstein condensates. Our 
results provide direct verification for the experimental 
observations of Bretin et al. 14[. We interpreted the 



coupling between the external perturbation and the exci- 
tation of Kelvin waves in terms of an elliptical instability 
mechanism driven by a parametric resonance between the 
counter-rotating quadrupole surfon and a pair of kelvons 
of opposite momenta. Subsequently, the Kelvin waves 
were found to decay to longer wavelength excitations via 
phonon emission. We presented a modification to the 
Kelvin wave dispersion relation for inhomogeneous su- 
perfluids. Finally, we proposed an experimental method 
to excite Kelvin waves of any wave number. We have 
shown that the pure kelvon decay provides a powerful 
damping mechanism in superfluid systems even at zero 
temperature. This may be viewed as an elementary de- 
cay channel for superfluid turbulence, the higher order 
nonlinear processes involving vortex (self) reconnections. 



Our methods also enable us to go beyond the in-plane 
Tkachenko modes 30|, 31, 32, 33, 34 1 and to study the full 
dispersion relation of Kelvin-Tkachenko collective modes 
in three-dimensional vortex lattices. 
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